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■^T . Abstract. We consider the problem of reconstructing a potential V 

from the scattering amplitude at a fixed positive energy. More specif- 
ically, we consider the method of Bukgheim to recover V from the 
Dirichlet-to-Neuman map on the boundary of a bounded planar domain 
which contains the support of the potential. We draw a connection 
between this and a question of Carleson regarding the convergence to 
initial data of solutions to time-dependent Schrodinger equations. As 
C^ ' a consequence we are able to recover the potentials with only half a 

"tl2 . derivative in L 2 . We also show that our result is sharp. 
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1. Introduction 



We consider the Schrodinger equation An = Vu on a bounded domain Q 
in the plane. For each solution u, we are given the value of both u and Vu • n 
on the boundary d£l, where n is the exterior unit normal on dSl. The goal 
is then to recover the potential V from this information. 

We suppose throughout that V £ L 2 (f2) and that is not a Dirichlet 
eigenvalue for the Hamiltonian — A + V. Then for each / € i7 1 ' 2 (9n), there 
is a unique solution u in H 1 ^) to the Dirichlet problem 



% ! (1) i-i« - /■ 



and the Dirichlet-to-Neumann (DN) map Ay can be formally defined by 

Ay : / i-+ Vu-n\QQ. 

Then a restatement of our goal is to recover V from knowledge of Ay. 

We come to this problem via a question of Calderon regarding impedance 
tomography [12] , however the DN map Ay_ K 2 and the scattering amplitude 
at energy n 2 are uniquely determined by each other, and indeed the DN 
map can be recovered from the scattering amplitude (see the appendix for 
explicit formulae) . Thus we are also addressing the question of whether it is 
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possible to recover the potential from the scattering data at a fixed positive 
energy. 

In higher dimensions, Sylvester and Uhlmann proved that smooth poten- 
tials are uniquely determined by the DN map [50] (see [37l [38l E] for non- 
smooth potentials and [USUI [22] for the conductivity problem). The unique- 
ness result was extended to a reconstruction procedure by Nachman [32l [33] . 
The planar case is quite different mathematically as it is not overdetermined. 
Here the first uniqueness and reconstruction algorithm was proved by Nach- 
man [34] via 9-methods for potentials of conductivity type (see also [10] for 
uniqueness with less regularity). Sun and Uhlmann \46\ 148] proved unique- 
ness for potentials satisfying nearness conditions to each other. Isakov and 
Nachman [25j then reconstructed the real valued L p -potentials, p > 1, in the 
case that their eigenvalues are strictly positive. The 3-method in combina- 
tion with the theory of quasiconformal maps gave the uniqueness result for 
the conductivity equation with measurable coefficients [2]. The problem for 
the general Schrodinger equation was solved only in 2008 by Bukhgeim [TTJ 
for C 1 -potentials. Bukhgeim's result has since been improved and extended 
to treat related inverse problems (see for example [391 I2H [23] ) ■ 

The aim of this article is to emphasise a surprising connection between 
the pioneering work of Bukhgeim [TTJ and Carleson's question [13] regarding 
the convergence to initial data of solutions to time-dependent Schrodinger 
equations. Elaborating on this new point of view we obtain a reconstruction 
theorem for general planar potentials with only half a derivative in L 2 , which 
is sharp with respect to the regularity. The precise statements are given in 
the forthcoming Corollary 11.31 and Theorem 11.41 

To describe the results in more detail, we recall that the starting point 
in [TTJ was to consider solutions to Am = Vu of the form u = e 1 ^ (l + w), 
where from now on 

ip(z) = ipk,x(z) = g{z — x) , z € C, x € Q. 

Solutions of this type have a long history (see for example [2TJ [501 123 HE] ) , 
and in this form they were considered first by Bukhgeim. We will recover 
the potential with a countable number of measurements on the boundary, 
so we take k £ N. We will require the homogeneous Sobolev spaces with 
norm given by H/H^p = ||(— A) s ' 2 /||^2, where (— A) s > 2 is defined via the 
Fourier transform as usual. In Section [21 we prove that if the zero extension 
of the potential V is contained in H s with < s < 1, and k is sufficiently 
large, then we can take w = Wk )X € H s with a bound for the norms which 
is decreasing to zero in k. We write Uk )X = e*^(l + w) for these w G H s . 

The definition of the DN map yields the basic integral formula in inverse 
problems; Alessandrini's identity . Indeed, if u, v E H 1 ^) satisfy Au = Vu 
and Av = 0, then the formula states that 

((A v - A )[u],v) := / (A v -A )[u]v= / Vuv. 
x ' Jan Jn 
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Taking u = Uk, x , which is also in H 1 (Q), and v = e 1 ^ this yields 

(2) ((A y - A )[u M ], e^ ) = J e^ + ^V(l + w) , 

and so the integral over f2 can be obtained from information on the boundary. 
The bulk of the article is concerned with recovering the potential from 
the integral on the right-hand side of ([2]). However, in order to calculate 
the value of the integral, without knowing the value of the potential V 
inside fi, we need to calculate the value of the left-hand side of ([2]). That 
is to say, we must determine the values of Uk, x on the boundary from the 
DN map. In the case of linear phase, this was achieved by Nachman [M] 
for L p -potentials V, with p > 1, and Lipschitz boundary. For C 1 -potentials, 
with C 2 -boundary, the result was extended by Novikov and Santacesaria 
to quadratic phases [55] . Here we show that for quadratic phases almost 
no regularity is needed. We consider potentials in the inhomogeneous L 2 - 
Sobolev space H s , defined as before with (— A) s ' 2 replaced by (I — A) s ' 2 . 
Our starting point is similar to [33] but we give a shorter argument, avoiding 
single layer potentials. 

Theorem 1.1. Let V G H s with s > and suppose that f2 is Lipschitz. 
Then we can identify bounded operators Rk '■ H 1 ' 2 ^^!) — > H ' 2 (d£l), de- 
pending only on k and Ay — Aq, such that 



Uk,x\dn = Rk[e 



iip 



dn\ 



For C 1 -potentials, Bukhgeim [TT] proved that the right-hand side of ([2]), 
multiplied by (4Tri)~ 1 k, converges to V{x) for all x € O, when k tends to 
infinity. In Section HI we obtain this convergence for potentials in H s with 
s > 1. For discontinuous potentials we are no longer able to recover at each 
point. Instead we bound the fractal dimension of the sets where the recovery 
fails. As Sobolev spaces are only defined modulo sets of zero Lebesgue 
measure, we consider first the potential spaces L s ' 2 = (— A)~ S ' 2 L 2 (M. 2 ), and 
bound the Hausdorff dimension of the points where the recovery fails. 

Theorem 1.2. Let V G L s ' 2 with 1/2 < s < 1. Then 
dim^jxe Q : ^( (Ay - A )[u M ],e^ -frV{x) as k -> oo} ^ 2 - s. 

As the members of H s coincide almost everywhere with members of L s ' 2 , 
we see that rough and unbounded potentials can be recovered almost every- 
where from information on the boundary. Note that these results are stable 
in the sense that k G N can be replaced by any sequence {nk}k€N such that 
rik tends to infinity as k tends to infinity. 

Corollary 1.3. Let V G hi 1 ! 2 . Then 

lim ^j/(A v - Ao)[u k ,x},e^) = V(x), a.e. x£Q. 
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In Section El we will prove that this is sharp in the sense of the following 
theorem. Note that even though there is divergence on a set of full Hausdorff 
dimension when s < 1/2, the dimension of the divergence set is bounded 
above by 3/2 when s ^ 1/2. 

Theorem 1.4. Let s < 1/2. Then there exists a potential V £ H s , sup- 
ported in Q, for which 

'xGO:^((Ay-Ao)[u M ],e ! ^>F(x) as fc^oo}|^0. 

Blasten [7] proved that potentials in H s with s > are uniquely de- 
termined by the DN map. It is a curious phenomenon that, within the 
Bukhgeim approach, uniqueness and reconstruction have different smooth- 
ness barriers. 

The DN map Ay_ K 2 can be recovered from the scattering amplitude at a 
fixed energy k 2 > (see the appendix), from which we are able to recover 
the potential V — n 2 Xtt rather than V. We are free to choose the domain ft. 
Taking 0, to be a square, we obtain the following recovery formula. Here 
Uk t x are Bukhgeim solutions which solve Ait = (V — k?)u in Q. 

Theorem 1.5. Let V € H 1 ' 2 be supported in a square O. Then 

lim ^-((A v _ k2 -A )[U kiX ],e^) + K 2 = V(x), a.e. x G fi. 

Interpreting the problem acoustically, it is unsurprising that we are unable 
to recover potentials in H s with s < 1/2. Taking 

V{x) = k 2 (1-c- 2 (x)), 

where c(x) denotes the speed of sound at x, the scattering solutions it also 
satisfy c 2 Au + k 2 u = 0. Now there are potentials in H s , with s < 1/2, which 
are singular on closed curves (see for example [52J ) . Thus the speed of sound 
is zero on the curve and so a continuous solution u would be zero. That is 
to say, the continuous incident waves cannot pass through the curve, and 
we should not expect to be able to detect a modification of the potential in 
the interior. From this point of view, the uniqueness result of Blasten [7] 
reflects tunneling phenomena in quantum mechanics. 

2. The Bukhgeim solutions 

Writing d z = \{d x — id y ) and c\ = \{d x + id y ), we consider the complex 
analytic interpretation of the Schrodinger equation Ad z c\u = Vu. When 
looking for solutions of the form it = e l ^(l + w) , the equation is equivalent 
to the system 

2c\w = e-^+^v, 2d z v = eW + ^V(l + w), 
which is solved in whenever 



lQ-1 



iv - jd- 



e-'^+^XQ d~ x [ e ^ + ^ V(l + w) 
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Here, we take Q to be a fixed, auxiliary, axis-parallel square which properly 
contains 0. Thus, defining the operator Sy = Sy X by 

S V [F] = Id- 1 le^+^XQ d- 1 [eW*® VF] 

we see that as soon as ||Sy-||jj s jp < 1, we can treat (I — Sy) -1 by Neumann 
series to deduce that it is a bounded operator on H s . This yields a solution 
Uk, x = e J ^(l + w) where 

(3) w = w k = Q.-Sk)- 1 S$[l]eH a . 

In what remains of this section, we prove that Sy is contractive for suffi- 
ciently large k. This property will be crucial in the proof of Theorem II .11 as 
well as in Section H We write S v [f] = \Si[Vf], where 

S k = d= l o M- k o d- 1 o M fc 

and the multiplier operators M ±k are defined by M ±fc [F] = e ±z ^ + ^XQF- 
The key ingredient in the proof of the following estimate, is the classical 
lemma of van der Corput 



Lemma 2.1. Let < A < s\, S2 < 1. Then 

l|M ±fc [F](-,x)||^_. a < Ck- x \\F(;x)\\^ sl , xen. 
Proof. We will first prove that 

(4) ||M ±fe F|U-i ^ Ck~ l \\F\\ ox , 

V ' " Mi( 2,oo " "-"2,1 

where the Besov norms are defined as usual by 

Here, P~? = t?(2^'| ■ |)/ with •& satisfying supptf C (1/2,2) and 

As ||F|| o-i ^ C'llFHoo and \\F ||i ^ C||F|| ai , the estimate flU) would follow 

2,oo 2,1 

from 

(5) ||M^>|U < Ck~ x \\F\\ X . 

Now, by the Fourier inversion formula and Fubini's theorem, 

|M±fc>(£)| = _!_ f e ±W(*)+tfW) ( ' F(u)e iz - w due- iz <dz 
( 27 Jq J 

so that ([5]) follows by two applications of van der Corput's lemma [18] (fac- 
torising the integral into the product of two integrals). 
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Now, by a trivial decomposition, for j = 1, 2, it will suffice to prove 
(6) IIMj^llH-a < Ck- X \\F\\^ S1 , < A < si, s 2 < 1, 

whenever F is frequency supported in D or M 2 \ID), where 

Mf k F = (m^Fxo) and M± k F = 

so we have four cases. It is clear that ([5]) and thus ([4]) also hold for these 
slightly modified operators. 

By summing geometric series in (j4]), we obtain 

with S4 less than or greater than one depending on whether we are in the 
case j = 1 or 2, and similarly for S3. On the other hand, by the Hardy- 
Littlewood-Sobolev inequality, Holder's inequality and interpolation, one 
has 

l|Mf[F]||^_ S6 <C||F||^ BI 0<*5,86<1. 
Interpolating between the two estimates yields ([6]), and so we are done. □ 

In the following lemma, we optimise the decay in k, which will be impor- 
tant in Section HI 

Lemma 2.2. Let < s,X < 1. Then 

\\S k [F](;x)\\ tis ^ Ck- X \\F(;x)\\ 6a , XeQ. 

Proof. By two applications of Lemma I2.1( 






^Ck- Xl \\M k \ 



where < Ai < 1 — s and < A2 < s, and so we are done. □ 

In the following lemma, we use Lemma 12.11 only once, and gain some 
integrability using the Hardy-Littlewood-Sobolev theorem. By taking k 
sufficiently large, we obtain our contraction and thus our Bukhgeim solution 
u = u k,x as described above. 

Lemma 2.3. Let < A < 2s, 1 - s. Then 

\\Sv[F](;x)\\hs < Ck- x \\V\\ A ,\\F(;x)\\z., xe!l. 
Proof. By the Cauchy-Schwarz and Hardy-Littlewood-Sobolev inequalities, 

\\VF\\ q ^ \\Vh q \\F\\ 2q ^C\\V\\^\\F\\^, 
where q = j^—. Thus, as Syfi 7 ] = S^[VF], it will suffice to prove that 

im fc ll ■ < Ck~^ 
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When < s < 1/3, by Lemma 12. 11 we have 

^C7A;- A ||M fc || L ^ L9 , 
where < A < 2s, 1 — s. When s ^ 1/3, we also use Holder's inequality; 
HWis-tH- < llM^oa^oM^II^^, ^C7A;- A ||a z - 1 oM fe || L ,^ 1 _ i 

^ OK ll iV1 IIl<3^_h"- s 

5% OK ||1V1 ||x,«_ + x,«* , 

where g* = ^rj, and so we are done. □ 

3. Proof of Theorem 11.11 

In this section we show that the boundary values of our Bukhgeim so- 
lution Uk x can be determined from knowledge of Ay. The argument is 
inspired by [34^ Theorem 5] but we replace the Faddeev green function Gk 
by its analogue in terms of the operator S v and avoid the use of single layer 
potentials. 

Indeed, considering the kernel representation of S^, we can write Sy[-F] 
in the form 

S v [F}(z)= f g^z,v)V(v)F{r,)dv. 

Jn 

where cju, the kernel of Sj, is given by 



m 
e 



g ^ z ,rj) = XQ (r ] r-—- 2 / — — ) e -(*(«)-HK«)) fa. 

47r JQ(u-r])(z-u}) 



In order to work directly with exponentially growing solutions we conju- 
gate g^h with the exponential factors, so that 



(7) / G^z,r,)V(r,)F(r,)d V = ^^S v [e^F](z), 
Jn 

where G^,(z,rj) = e 1 ^^ g^(z,r])e~ t ^ v \ Notice also that when z G Q\£l and 
7] G f2, we have that 

A v G i ,(z,r)) = 0. 
Thus, if we take ((7J) with F = Py(f), where Pv(f) solves Au = Vu with 
u\dn = f, using Alessandrini's identity we obtain that, for each z G Q \ O, 

(8) ((Ay - Ao)[f],G^(z,.)\ dn ) = e^Syie-^PviMz). 

In particular the right-hand side belongs to H 1 (Q \ O) and hence we can 
define the operator Ew, : H 1 ' 2 — > H 1 ' 2 by 

r^[f]=T r o((A v -A )[f},G t \ dCl 
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where T r : -ff 1 (Q \ f2) — » i?" 1 / 2 (3fi) is the trace operator. Now, by the 
definitions of Uk, x and w, we also deduce from and (|3|) that 

(9) [ G^,ri)V{fi) UkM dv = e**S$r[l + w] = e^w = u M - e**. 
Combining fl7|, ([5]). and @ we obtain the integral identity 

(I - r^)[«fe )a .|an] = e \aci- 

Thus, we can determine Uk X ° n the boundary if we can invert (I — 1^,). 
By the Fredholm alternative it will suffice to show that 1^, is compact and 
that (I — FU) has a trivial kernel on H 1 ' 2 (dQ). 

Theorem 3.1. Let V € H s with < s < 1. Then 
(i) r^ is compact 

(ii) l>[/] = /^/ = o. 

Proof of (i). We have that 

T4f]=T r [e^S k v [e-^P v (f)}]. 

As the set of compact operators is a left and right ideal, we consider the 
boundedness properties of each component of the composition. Firstly, Py : 
H^ 2 (dQ) -»■ H x {n) is bounded. Secondly, H 1 ^) ^ IP (SI) compactly for 
all 2 < p < oo. Now taking p sufficiently large and ^ = — I — > by the 
boundedness of the Cauchy transform followed by the Hardy-Littlewood- 
Sobolev inequality, 

\\Svi e G]\\ H i( Q \ty ^ C\\VG\\ L 2^ < C||y|| L? (n)||G'||LP(Q) 

^C\\V\\ H3 \\G\\ LP{Q) . 

Finally, T r : H X (Q \ U) ->• F 1 / 2 ^) is bounded. Since the embedding 
if 1 (0) ^-t- L P (Q) is compact, it follows that Ew, is compact. 

Proo/ o/ (mJ. Letting p = Sy[e~^Py(/)], we have that 

c\[e^p] = le^XQd-^VPvif)], 

so that 

4d x d g [e i *p] = VPv(f) on fi. 

This can be rewritten as A[e 2 ^/o — Py(/)] = on O. Now by hypothesis 
I^,[/] = /, so that by (|5J) we have e 1 ^ p = / on 90. Combining the two, we 
see that 

e ^p = p v (f) on ft. 

From the definition of p we see that p = S v [p], and as soon as S v is strictly 
contractive, that p = 0. This of course follows from Lemma 12.31 for large 
enough k. Thus, / = e 1 ^ p = 0, so that I — L^, is injective as desired. □ 
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4. Potential recovery 

In order to recover the potential at x € CI, it remains to show that the 
right-hand side of Alessandrini's identity ([2]) converges to V{x). That is to 
say T\ +w V(x) converges to V(x) as k tends to infinity, where 

T\ +w [F](x) = A f e W«+^> F(z) (i + w{z) ) dz . 

First we show that T^V can be considered to be a remainder term. 
Theorem 4.1. Let V € H s with < s < 1. JTien 

lim T*[V](x)=0, ie!l. 

fe— >oo 

Moreover, i/ fe ^ (1 + c\\V\\ 2 ^ 3 ) max{ ^' T ^ } and < A < s, i/ien 

sup|T^y](x)KCA;- A ||y||^. 

Proof. By Lemma 12. 11 

iT^Vl^l^CfcllM^II^IMI^ 

where < A < s. By Lemma l2.3| we can treat (I — S^) _1 by Neumann 
series to deduce that it is a bounded operator on H s . In particular, when 
s > 1/3, this is true if k ^ 1 and Ck ~ ||V[|#« < \. Then 

\^t\V\{x)\^C^- x \\V\\ tis \\^\\V\\\ tia 

^Ck- x+e \\V\\% s , e>0, 

by an application of Lemma 12.21 which is the desired estimate. □ 

Noting that e i(^)+^)) = eX p (ik {zi ~ Xl)2 ~^ Z2 ~ X2)2 ), it remains to prove 
(10) ]imT ! l[V](x) = V(x), 

where T^ is defined by 

Jin 



T*[F](x) = -f- I exp ( ifc (*i-*i)M*2-* 2 ) 2 ) F ( Z )^, 
4iri J 

Now when i 7 is a Schwartz function, this is equal to e** F(x), where 

(2tt) 2 7 M 2 
This follows easily, making use of the distributional formula 
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which holds for Schwartz functions (p. We see that when V is a Schwartz 
function, T\V solves the time-dependent nonelliptic Schrodinger equation, 

idtu + Du = 0, 

where □ = d xlXl — d X2X2 , with initial data V at time 1/k. When V G H s 
with s > 1, both V and its Fourier transform are integrable, and so both 
T\V and e*fc V are continuous functions which are again equal pointwise. 
Thus, in the following lemma we obtain the convergence (I10p and therefore 
complete the reconstruction for potentials in H s with s > 1. 



Lemma 4.2. 



Let V E H s 



with 1 < s < 3. Then 
*V(x)-V(x)\ <Ck^\\V\\ H ; x € fi. 
Proof. By the Fourier inversion formula and the Cauchy-Schwarz inequality, 



le* 



„itn 



V(a;)-V(ar)| 



1 



$; 



(2vr) 2 

ll^llfl 



*W P -*ir«?-3) 



V(£)e'"*'~(e "« 



1 



IVI 



2fc" 



i? s 



l£l 2 * 
2-2cos(i(^ 






1/2 



Kl 



2* 



^ 



1/2 



< 2AT 



1^1 



1^1 



ff s 



ff s 



ain 2 Q(tf-& 2 )) sV 2 



|£|2(s-2) 



|£|b 



^ 



1 \l/2 



2..s 



where we have used the trigonometric identity 2 sin 9 
fact that sin# ^ |0|. 



Altogether we see that |T* 



cos 20 and the 
□ 



l+w V(a>) - V(a?)| < Ck~ for all i € and 
V € -£P with 1 < s < 3, which improves upon the decay rate of [39] where 
they recovered C 2 potentials. Note that there can be no decay rates, at least 
for the main term, for the potentials of H s with s ^ 1 as they would then 
be uniform limits of continuous functions and thus continuous. 

For discontinuous potentials we are no longer able to recover at each point. 
Instead we bound the fractal dimension of the sets where the recovery fails. 
This point of view has its origins in the work of Beurling who bounded 
the capacity of the divergence sets of Fourier series [6] (see also [3]). Now 
Sobolev spaces are only defined modulo sets of zero Lebesgue measure, and 
so we consider first the potential spaces 



L s > z = {I s *g : g£L 2 (I 

where I s is the Riesz potential | • | s_2 . As I s (£) = 
is also a member of (an equivalence class of) H 



CM 



5 , we have that I s * g 
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To bound the dimension of the sets where the recovery fails, we will prove 
maximal estimates with respect to fractal measures. We say that a positive 
Borel measure \x is a-dimensional if 

(11) ca(ji) ■= sup ^ V V " < oo, < a < 2, 

x£K 2 ,r>0 r 

and denote by A4 a (£l) the a-dimensional probability measures which are 
supported in $7. For < s < 1, we will require the elementary inequality 



( 12 ) Us * gW&idn) ~ VCoG") lbllz,2 (R 2), a > 2 - 2s, 

which holds whenever /x G _M a (f2) and g G L 2 (M 2 ). To see this, we note 
that by Fubini's theorem and the Cauchy-Schwarz inequality, 

\\h *g\\L i (dn) < ll-Ts *HIl 2 IIs , IIl 2 , 

so that (fT2l) follows by proving 

HJs *Hli 2 < Ca(ju), a>2-2s. 
Now by Plancherel's theorem, 

\\i s * t i\\l 2 = (2n)-^T s Ml 2 <J^omh s (0^<J^*h s (.y)d^y) 

/" /" dfj,(x)dfi(y) 
J J \x-y\ 2 ~ 2s ' 

which is nothing more than the (2 — 2s)-energy. Then, by an appropriate 
dyadic decomposition, 



whenever a > 2 — 2s and fi G Al a (0). 

The Fourier transform of less regular potentials V is not necessarily inte- 
grable, and so in that case e J fc V is not even well-defined. Instead we make 
do with the pointwise limit 

(13) Tf M(x) = lim G N * T$[V](x) = lim e^ n [G N *V](x), x£ ft, 

N— s>oo " ' N— s>oo 

where Gn = N 2 G(N-) and G is the Gaussian e~'' . This formula holds as V 
is compactly supported and integrable; conditions which the initial data in 
the time-dependent theory does not normally satisfy. We will also require 
the following lemma due, in this form, to Sjolin [43J. 

Lemma 4.3. [43] Let x, t G R, 7 G [1/2, 1) and N ^ 1. Then 



<^_ 



n{N- l £) e*^-* 8 ) ,. 

KF * 

where the constant implied by the symbol < depends only on 7 and the 
Schwartz function 77. 
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In the following theorem, we employ the Kolmogorov-Seliverstov-Plessner 
method, as used by Carleson [13] for the one-dimensional Schrodinger equa- 
tion. Dahlberg and Kenig [19] proved that the result of Carleson is sharp and 
noted that his argument could be applied to the higher dimensional problem 
(for which the argument is no longer sharp for the elliptic equation, see [8]). 
We refine their argument, which extends to the nonelliptic case, by proving 
estimates which hold uniformly with respect to fractal measures. 



Theorem 4.4. Let 1/2 ^ s < 1. Then 
^ n [G N *I s *g]\\\ TXMiA < 

k^l JV>1 

whenever fi £ A4 a (Q) and g G L? . 



swpsup\e l k u [G N *I s *g]\\\ Ll{dfi) < v/c« (a*) IMI */»(*»), 



a > 2 



Proof. By linearising, it will suffice to prove 



(14) 



Jt(w)D, 



< '" ~[Gn(x) *I s *g] w(x) d/i(x) 



< c^C/j,) Il^llia, a>2 



uniformly in measurable functions £ : f2 — >• R, iV : SI — >■ N and w : CI — > B. 
By Fubini's theorem and the Cauchy-Schwarz inequality, the left-hand side 
of (fl4"l) is bounded by 

2 



l?(OIX 



G 



£ 



JV(x) 



M*m-&) e ix <w(x)dft(x) 



_d£_ 



Writing the squared integral as a double integral, and applying Fubini's 
theorem again, it will suffice to show that 



G 



Z \ „( i \ e i(t(x)~t(y))(e i -e 2 ) e i(x-y)<_ d Z 



N(x 



G 



N{x) 



iei 2s 

(15) w(x)w(y) d/j,(x)d/j,(y) < c a (fj,) 

uniformly in the functions t, N and w. Now, as |£| 2s ^ |£i| s |£2| S j the left- 
hand side of (fT5|) is bounded by 

2 



n 



G 



i, 



G 



i 



., \ c i{-VP+ 1 Wx)-t[y))e ic i{x i -y i )i 1 d ^3 



-N{x)J \N(y)J fo| 

w(x)w(y) dfi{x)dfi{y), 
and by Lemma 14.31 we have 



f 



101 



-d& 



< 

no 



y.) 



.11— s 



Substituting in, we see that the left-hand side of (|15p is bounded by 



(16) C 



\w(x)w(y)\dfi(x)dfi(y) 
\%l -yi| x ~ s k2 — J/2| 1 ~ s 



^C 



d/j,(x)dfj,(y) 



\xi -yil 1 s \x 2 -2/2I 1 s 
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To complete the proof, we are required to bound (TT6j) by c a (fi). This will 
require a dyadic decomposition which lends itself to the singularities along 
the axis-parallel lines A y defined by 

A y = {x G Q : xi = y\ or X2 = yi\ , y G £1- 

Covering A y by balls {Bj}j^i of radius r,j and using the definition (fTTj) of 
c a (/j,), we have 

Taking the infimum over all such coverings and using the fact that the a- 
Hausdorff measure of A y is zero when a > 1, we see that (J,(A y ) = for all 
fj, G .M a (f2). Thus we can ignore the sets A y when decomposing the inner 
integral of (fT6|) . 

For each j, £ G Z we break up Q D fi into dyadic rectangles of dimensions 
2~ J x 2 and consider the unique rectangle Rj/ which contains y. We call 
the unique rectangles Rj-\^-\, Rj-i t £, and Rj,£-i that contain Rj t £, the 
mother, the father, and the stepfather respectively. We write Kfg ~ Rjj if 
their mothers touch, but their fathers and stepfathers do not. As /J,(A y ) = 0, 
we can write 






F{x,y)dn{x), 



i,O0 n:R" ~JL, " R j,l 



















































































■y 



















































































The rectangles of dimensions 2 J X 2 , with 1 ^ j,£. ^ 3, 
associated with a single point y. 



which yields 



Without loss of generality, we can suppose that 

£>j^0 n:R n e ~R jt e j^>0 n:R n ^R^/, 
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so that 

® < c f E E ^-'W-^mdM- 

Now by covering each rectangle by discs of radius 2~ 3 , and using the defini- 
tion (|11|) of c a (p) , we see that 

and for each rectangle Rj t £ there are exactly nine rectangles R% which satisfy 
R ],e~ R j,e- Thus 

(USD < cM £ 2* 2 — a >2-* < cM, 

when a > 2 — s, and so we are done. □ 

Proof of Theorem ] 1.2k By Alessandrini's identity (|2j) and Frostman's lemma 
(see for example [3XJ ) , it will suffice to prove that 

(17) Jx : limsup \T k l+w [V](x) - V{x)\ / o} = 

*• fc— >oo ' 

whenever \x G M a (Q) and V G L s ' 2 (fi) with a > 2 - s. By Theorem |4"TT1 
and (|13p . this would follow from 

/ijx : limsup limsup |e** D [Grjv * V](x) - V(x)\ ^ o| = 0. 

Writing V = I s * g, where g G L 2 , we take a Schwartz function /i so that 
\\g — h\\i'2 < e. Then 



fil x : limsup limsup \e l k U [G N * V](x) - V(x)\ > A \ 

*• k— >oo N^t-oo ^ 

^ fil x : sup sup |e l * p [Gjv * I s * (g - h)](x)\ > A/3 } + 
*- fe>i jv>i - 1 

//^ a; : limsup limsup |e l fc D [G/v * 7 S * h](x) — I s * h(x)\ > A/3 > + 

*■ k— >-oo TV— »oo ' 

n{x : |I S *(/*-<?) (x)| >A/3}. 

As the terms involving /i are continuous in all parameters, the second set of 
the three is empty, so by the elementary inequality (|12p and Theorem 14,41 
we see that 

fj,\x : limsup \T\ +w V(x) -V(x)\ > \\ < A -1 y/cJjJL) \\g - h\\ L 2 

< A"V c a(^)e, 
for all e > 0, which yields (|17fl . and so we are done. □ 

Proof of Theorem \1.5l For V G if 1 ' 2 , the potentials q = V — K 2 xn are 
contained in H s for < s < 1/2 (see for example [20J) and so we find 
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Bukhgiem solutions Uk, x , associated to q, and recover their value on the 
boundary as before. Corollary 11.31 requires the potential q to be contained 
in H 1 ' 2 which is not satisfied for any domain. However, it is clear from the 
proof of Theorem 14.41 that we can relax this additional condition further to 

IK^) 1/4 (^ 1/4(Z II^(K 2 ) <00 ' 

which is satisfied when fHs a square, but not when it is a disc. □ 

Finally we remark that the uniqueness result of Blasten [7j can be ob- 
served using the connection with the time-dependent Schrodinger equation. 
Indeed if the scattering data or boundary measurements are the same for 
two potentials V\ and V 2 , then by Alessandrini's identity (J2J), 

\\V 2 - V x \\ L i = \\V 2 - T k 1+w V 2 + T k x+w V x - ViHra, 

so that by the triangle inequality and Lemma 14. 1\ it suffices to prove 

\\V -T^V\\ L 2 ^0 as k ^00, 

which is a well-known property of the Schrodinger flow. 

5. Proof of Theorem 11.41 

First we construct a real potential V, supported in $7, and contained in 
H s with s < 1/2, for which 

(ie!l: lim e^ D [V]{x) ■/> V(x)\ + 0. 

I fc— >oo J 

Throughout this section we work with a different set of coordinates from the 
previous sections. Indeed, for Schwartz functions F, we now write 

el4D[i?](x) = Wfj eix<e ~ m ^ 2 %) <%• 

Let (j) be a positive Schwartz function, compactly supported in [—1/4, 1/4], 
and consider cfi = (f> * ip , which is supported in [—1/2,1/2]. Note that 
4> = (^o) 2 ^ 0. We consider the potential V defined by 

V(x) = S £ J V j {x) = s £ j 2 i < l -^ +l cos{2ix2)(t>{2 j x 1 )<i){x2) 

= J^2( 1 -® j J 2JX2 <t>(2 j x 1 )<f>(x 2 ) + ^2( 1 - /3 ^e- i2 ^ 2 < /»(2 J xi)^(x 2 ) 

j>2 j>2 

j>2 j>2 

which is supported in [— |, |] x [— i, ^]. If /3 € (1/2 + s, 1), by changes of 
variables, 

\\v\\ 2 HS < c^2< 1 - 2 ^ /i^i)tei 2 (i+iei 2 m < oo. 



Hi 
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Thus V is finite almost everywhere, and we will show that e 1 * U V diverges 
°n[^,|]x[-^,^]. 

This potential is an adaptation of an initial datum for the time-dependent 
nonelliptic Schrodinger equation considered in [41] . The initial datum there 
was not real, the diverging sequence of time was allowed to depend on the 
point x, and more crucially, the initial datum was not compactly supported. 
Thus our arguments will have a different flavour, working on the frequency 
and spatial side simultaneously. 

By changes of variables and the Fourier inversion formula, 

e* in [U+](x) = 2(1 "f'f X2 / 0(6)0(6) e-^^V^i^-^'H^)^ 



(2tt) 2 
2?r 



4>{2\ Xl - y+H - 2t£ 2 ))0(£ 2 )e^ 2 d6- 



Taking t = l/k with k the nearest natural number to 2 J+ /x%, 

2 (l-P)J e i2Jx 2 



e^ u [V+](x) 



2tt 



HC(^J)- n^6)0(6) e** x *dZ2, 



where \({xi,j)\ ^ \ when x\ € [jq, j], so that, using the compact support 
of 0, we see that 



\e^ u [V+]{x)\ 



2 (l-j8)i /-16 



2tt 



-16 



0(C(*i,i)- ^6)0(6) e**°»dS2 



> 



2 (l-/3)j /•!<• 



2?r 



16 



0(C(»i,i) - n^6)0(6) cos(c 2 x 2 )d6 



Now when x 2 € [— -nr, -nr], we have |6 X 2| ^ 1> so that | cos(£ 2 x 2 )| > cos(l) 
Using the fact that and <p are nonnegative, we obtain 

2( 1 -Wcos(l) rW 

J -16 



\e^ D [V;](x)\ > 



2tt 



0(C(*i,i)- ^6)0(6)^6 



It remains to bound from above the solution associated to the other pieces 
of the potential. Again, by the Fourier inversion formula, 

|e^ D [V±](*)| " 



(2tt)2 
2 (i-/9)/ 



*+i. 



0(6)^(6) e^ 1 ^ 2 &(*it— J+k 3 *)^ 



2tt 



(2 £ (x lT ^-|6))0(6)e^ 2 (i6 



Using the fact that 0(y) ^C\y\ 1 ' 2 , we obtain 



^[^Wl^^ 2 ^ 



10(6)1 



\X1 qp 2f±l _ 2 t 11/2 



d6- 
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Taking < e < min{l/4, 1 — f3}, and using the rapid decay of cj>, we see that 

\e^[V±](x)\ ^ C2(V»-W( f } dt 2 + C2-i 

W|&|<2«* |xiT^-|6| 1/2 

Now one can check that when £ 7^ j or j = £ and =p is an addition, 

1 9 j. q 9^+1 1 

l|6l^|kiTVl 

when |^2 1 ^ 2 je . Indeed, when j > £, the left-hand side is less than j\x\\ 
which is less than the right-hand side. On the other hand, when j < £ or 
j = £ and =F is an addition, the left-hand side is less than i|xi| which is 
less than the right-hand side. Thus, the integrand of the final integral is 
nonsingular so that the integral is bounded by C\xi\~ l ' 2 2 e 'J ^ C2 ej . 
By summing a geometric series in £, we obtain 

< C 2 2 ej , 



5>«i 


3 V e ± (x) + e t 


3 


» 


m 










le that 


on 


16' l] X 


L 16' 


16 J' 



|e^ D [y]| ^ |e^ D [^+]| - I JV* ^] + e^ D yr(x)| > Ci2^ 1 ~ /3 ) - C 2 2 je , 

which diverges as j tends to infinity. Considering forty-five degree rotations 
of the Vj, which are Schwartz functions, via the pointwise equality, this 
yields 

\T\[V}\ > |T£[K+]| - I ^Tj[Vf] +Tf[^"]| > dVQ-® - C 2 2 jt 

on a forty-five degree rotation of [^?, |] x [— i, ^?], so that |T*[V]| diverges 
as fc tends to infinity. Thus, by Theorem 14.11 combined with Alessandrini's 
identity ©, 

|x : s?/(Av-A )[u fc>a .|8n],e^|an)7^ 1 ^(a?) as fc -> °°} 

contains a forty-five degree rotation of [^, |] x [— y-r, i], which has nonzero 
Lebesgue measure. □ 

Note that this result is stable in the sense that k G N can be replaced by 
any sequence {nfc}fc g N satisfying n^ G [A;, A; + 1). 



Remark 5.1. In [44J, Sjolin asked for which values of s is it true that 
lim e^ A f(x) = 0, a.e. x G M d \(supp/), 

for all / G -SP. In principle, this question could have stronger positive results 
and weaker negative results than Carleson's question: for which values of s 
is it true that 

lim e^ A /(x) = f(x), a.e.xeM. d , 

k—>oo 
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for all / G H s ? Indeed, before Bourgain's recent breakthrough [8j, Sjolin 
proved a stronger positive result for his question than what was known for 
Carleson's question in three dimensions. Here we solve Sjolin's question 
completely for the nonelliptic equation in two dimensions. That is to say, 

lim e l k n f(x) = 0, a.e. x G M 2 \(supp/), 
for all / G H s if and only if s ^ 1/2. 

Appendix A. The DN map from the scattering amplitude 

It is well-known that in the absence of zero Dirichlet eigenvalues there is 
a unique weak solution to the Dirichlet problem (JTJ) that satisfies 

( 18 ) IMI#i(f>) ^ c \\f\\m/*{dQ.) 

(see for example |17| - in two dimensions L n ' 2 {W l ) can be replaced by L 2 (M 2 )). 
Here H l l 2 {dti) := if 1 (£})/#<} (ft), where i^(fi) denotes the closure of Cg°(fi) 
in H 1 ^). The DN map Ay is then defined by 

A v [f],i/j)= [ Ay[/]^= [vuV + Vu-VV, 

1 Jan Jn 

for all \E r G i? 1 (0) with ijj = ty + Hq(Q). When the solution and boundary 
are sufficiently regular, this definition coincides with that of the introduction 
by Green's formula. To see that Ay maps from i7 1 ' 2 (3Jl) to H~ 1 ' 2 (dQ), 
the dual of H l ' 2 {dSl), we note that by Holder's inequality and the Hardy- 
Littlewood-Sobolev inequality, 

A v[f],i>) ^ \\ u \\m(n)\\^\\m(n) + ll^ll2|k|| L 4 (Q) ||^|| L 4 (n) 
^ (I + C\\V\\ 2 )\\u\\ HHn) \\y\\ H i (n) 
whenever \I/ G iJ 1 (0), so that by (fTB"|) . we obtain 



Av[f],il>) ^C(l + \\V\\ 2 )\\f\\ Hl/Hdn) U\\ Hl/Hdn) . 

There are a number of different approaches to showing that the scatter- 
ing amplitude at a fixed energy k 2 > uniquely determines the DN map 
A V _ K 2 and vice versa (see for example [5j[32j[5TJH71ll9]). Here we follow a 
constructive argument due to Nachman [331 Section 3] . Instead of assuming 
that zero is not a Dirichlet eigenvalue of — A + V, we require now that k 2 is 
not an eigenvalue. This can be arranged by taking £1 sufficiently large as the 
eigenvalues decrease strictly as the domain grows [36J (the result of [30] can 
be extended to L 2 -potentials using the unique continuation of [26]). Here, 
unlike in the rest of the article, we additionally suppose that V is real. 

Let Gy and Go be the outgoing Green's functions that satisfy 

(-A + V- k 2 )G v (x, y) = 6{x - y), (-A - k 2 )G (x, y) = 6{x - y), 
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and let Sy and So be the corresponding near-field operators denned via 
single layer potentials; 



Sy[f](x) = G v (x,y)f(y)dy, S [f](x) = / G (x,y)f(y)dy. 
Jan Jon 

These are bounded and invertible, mapping H~ 1 ' 2 (d^l) to i? 1 ' 2 (9Q) (the 
two-dimensional proof can be found in [251 Proposition A.l]). Then Nach- 
man's formula [32], 

V-K 2 = A_ K 2 + Sy — S Q , 

allows us to recover the DN map. 

Thus it remains to recover the single layer potential Sy from the scattering 
amplitude Ay at energy k 2 . For weS 1 , the outgoing scattering solution 
v(-,uj,k) is the unique solution to the Lippmann-Schwinger equation 



(19) v(y,u,K) =e lKy -"- / G (y,z)V(z)v(z,u,K)dz. 

For (cr,(jj) € S 1 x S 1 , the scattering amplitude then satisfies 

(20) A v ((t,uj,k)= I e- iKa - z V(z)v{z,uj,K)dz. 



When Q is a disc, Nachman recovers Sy via formulae given by expansions 
in spherical harmonics as below. Otherwise he uses a density argument (we 
remark that Sylvester [39] also invokes density in order to recover). Since 
we have been obliged to work with $7 a square, at this point we deviate and 
instead follow an argument of Stefanov |45j . obtaining an explicit formula 
for the Green's function Gy in terms of Ay. Alternatively it seems likely 
that we could pass to the DN map on the square from that on the disc via 
the argument in [HU Section 6] for the conductivity problem, but we prefer 
this more direct approach. 

Stefanov worked in three dimensions, with bounded potentials, and a 
number of details change in two dimensions, so we present the argument. 
We recover Gy outside of a disc which contains the potential, but which is 
contained in the domain, so that Sy can be obtained by integrating along 
the sides of our square f2. First we require the following asymmptotics. 



Lemma A.l. 



-i e tK ' x ' e lK ' y ' ( x y 



Gy(x,y) - G Q (x,y) = j -Ay (— -,- — ,« 



o i a iJl~~ v 'VI I 7 I I 7 " " J ' ~ V i i — i i — 

" 7r/s: \x\ 2 \y\ 2 \ x \ \y\ \x p \y p 

Proof. It is well-known (see for example (3.66) in [35]) that Gy satisfies 

(21) G v (x,z) = - — -r v U,-7-r .« +° (— t ). 

(87r)2 K2 \x a v PI ' v \x a ' 
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and, in particular, 

(22) G (y, z) = G \ e-™m- z + o(-^ 

(87TJ2 K2 \y\2 \y\ 2 

From the definitions, it easy to verify that 

(23) G v (x, y) - G (x, y) = - [ G v (x, z)V(z)G (y, z) dz. 
Substituting in ([21]) and (j22j) . see that Gy{x,y) — Go(x,y) is equal to 

I r / e ^' z V(z)v (z, -|-| ,k )dz + o{ — i r , 



O/lft, £2 \y 2 J l x l X 2 \y 2 



so that (|20p and the symmetries of Ay yield the result. □ 

In the following, J n and Hn denote the Bessel and Hankel functions of 
the first kind of nth order, respectively (see for example |29j). We also write 
x in polar coordinates as (\x\,4> x ). 

Theorem A. 2. Let V € H s with s > be supported in the disc of radius p, 
centred at the origin, and consider the Fourier series 



A v (<t,u),k) = ^2 5Z 

nGZmGZ 



inipo- irrupt 



Then 



G v (x,y)-G Q (x,y) =£ ^-^e^a n , m H^( K \x\)H^( K \y\)e m ^e^y , 

where the series is uniformly, absolutely convergent for \x\ > \y\ > R > p. 
Proof. We can expand Go(x,y) = jHq (k\x — y\) as 

G (x,y) = H^\ K \x\)J ( K \y\) + 2Y,HJ l 1 \K\x\)J n ( K \y\)cos( ( l )x - (/),), 

(see for example [13 Section 3.4] or [U Theorem 3.4]). As H { } ] n = {-l) n H { n ] 
and J- n = (— l) n J n , in order to separate variables it will be convenient to 
write this as 

G (x,y) =Y / H0)(K\x\)J n (K\y\)e in ^e- in ^. 

neZ 

As before, it is easy to check that 

Gy(x, y) - G (x, y) = - G (x, z)V(z)G v (z, y) dz, 
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and so substituting (|23j) into this we obtain Gy — Gq = —I\ + I2, where 
h= [ G {x,z)V{z)G (z,v)dz 

h = G (x,zi)V(zi) G v (zi,Z2)V(z2)Go(y,Z2)dz 2 dzi. 



Now in both integrals we introduce the expansion of Go (note that Gq(x, y) = 
Go(y, x)), extracting the terms independent of z, Z\,z%. In this way we get 

(24) J i = E E ^ m H^(K\x\)H^(K\y\)e m ^e im ^, 

n£X m£Z 

(25) I 2 = E E ^n, m H^(K\x\)HW( K \y\)e m ^e im ^, 

where 
a 



V(z)J n (K\z\)J m ( K \z\)e-^ n+m ^ dz, 
P„,m = I Jn(K\z 1 \)V(z 1 )Gv(z 1 ,Z2)V(z 2 )J m (K\z2\)e- in ^e- im ^ dz x dz 2 . 



l n,m 



It remains to show that the sums (|24p and ()25[) converge uniformly and 
absolutely for |x| > \y\ > R > p. Once we know that this is is the case, we 
can take limits and use the asymptotics of the Hankel functions for large r; 

#(i) (r) = e -^i+4)(A)v + (J_ 



TIT 



7*2 



(see for example [29] Section 5.16]), and then Proposition I A . 1 1 tells us that 

2 ( — I) 171 



V f-" I li.l I ^ " I V.I I L I r\ — lb* lib 

IT 07T 

To see that the sums converge note that, by Holder's inequality, we have 

|"n,m| ^ Cpll^llL 2 ll^n( K l ' I) \\l°°(B p ) ||^m(«| ' \)\\l°°(B p )i 

\Pn,m\ < l|G'v||L 2 ( J B p x J B p )ll^ / llL2||«/n(K| • |) ||l°°(B p ) II ^m(«| ' |)||l,°°(J9 p )- 

At this point we deviate from [45J as there seems to be less local knowledge 
regarding Gy in two dimensions. Instead we can rewrite (|23|) as 

Gv(;y) = G (-,y) - (-A + V - k 2 - iO)" 1 [FG (-,y)], 

and use that the resolvent is bounded from L 2 ((l + | • \ 2 ) s ) to L 2 ((l + | • | 2 ) _<5 ) 
with 5 > 1/2 (see [TJ Theorem 4.2]). Thus, using that V is compactly 
supported, and taking ^ = ~ + g with large p so that 1 — - = s, 

\\Gv(;y)hHB p ) < \\Go(;y)\\mB p )+C p \\VG (;y)\\ L 2 {Bp) 

^\\G (;y)\\ LHBp) + C p \\V\\ q \\Go(;y)\\ LP{Bp ) 
<\\Go(;y)\\ L 2 {Bp) + C p \\V\\ H s\\Go(;y)\\ LP{ B p ), 
by the Hardy-Littlewood-Sobolev inequality. Integrating again with respect 
to y, and recalling that the singularity of Hq at the origin is logarithmic, 
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we see that ||GV||£2( B xB \ ^ C. Then, using the Taylor series expansion 
for the Bessel function, 



(-1)3 /r\2j+|n| 



I *-^ v.{\n\ + 7) V2 



lnl!V2y 



^Cprrrr(^) , < r < p, 



we see that 



|On,m| < C plrllL 2 TTr n 



n|!\2/ |m|!V2 

,, ,,q 1 1 0\\ n \ 1 / D\\ m \ 
Pn,m ^C p Vf Hs — (£) — (£) . 
|n|! V2/ |m|! \2/ 

Finally we require the Hankel function estimate, 

Ifl-pwKCflHid) 1 " 1 , R^r, 

which again follows from the power series expansion given that the function 
decays at infinity (see for example |291 Section 5.5] or [H Lemma 1]). The 
sums (|24l) and ([25]) are then bounded by a constant multiple of 

provided that |x| > \y\ > R > p, which is convergent, and so we are done. 

□ 
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